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^ ■ Abstract 

QQ , For non-zero I values, we present an analytical solution of the radial Schrodinger equation for the 

rotating Morse potential using the Pekeris approximation within the framework of the Asymptotic 

pg ■ Iteration Method. The bound state energy eigenvalues and corresponding wave functions are 

f^ . obtained for a number of diatomic molecules and the results are compared with the findings of the 

^^ super-symmetry, the hypervirial perturbation, the Nikiforov-Uvarov, the variational, the shifted 

\D 

^D ' 1/N and the modified shifted 1/N expansion methods. 
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analytical solution. 



I. INTRODUCTION 

The Morse potential has raised a great deal of interest over the years and has been one 
of the most useful models to describe the interaction between two atoms in a diatomic 
molecule. It is known that the radial Schrodinger equation for this potential can be solved 
exactly when the orbital angular quantum number I is equal to zero [1] . On the other hand, 
it is also known that for £ 7^ 0, one has to use some approximations to find analytical or 
semi-analytical solutions. Several schemes have been presented for obtaining approximate 



solutions 



Among these approximations, the most widely used and convenient one is the 



Pekeris approximation J3|,y], which is based on the expansion of the centrifugal barrier in a 
series of exponentials depending on the internuclear distance up to the second order. Other 
approximations have also been developed to find better analytical formulas for the rotating 
Morse potential. However, all these approximations other than the Pekeris one require 
the calculation of a state-dependent internuclear distance through the numerical solutions 
of .a„scende„ta. equation, fl B fl Q. I„ tK. .e.p.o, the .o.a«,, Mo.se pote,.,a. Kas 
SO far been solved by the super-symmetry (SUSY) [^, |9|, the Nikiforov-Uvarov method 
(NU) |lO|, |ll|, the shifted and modified shifted 1/N expansion methods jsj, |l2] as well as 
the variational method |l3| using Pekeris approximations for £ 7^ 0. It is also solved by 
using the hypervirial perturbation method (HV) [1J| with the full potential without Pekeris 
approximation. 

In this paper, our aim is to solve the rotating Morse potential using a different and more 



practical method called, the asymptotic iteration method (AIM) |l5l . Il6l | within the Pekeris 
approximation and to obtain the energy eigenvalues and corresponding eigenfunctions. In 
the next section, the asymptotic iteration method (AIM) is introduced. Then, in section UTTl 
the Schrodinger equation is solved by the asymptotic iteration method with the non-zero 
angular momentum quantum numbers for the rotating Morse potential: The exact energy 
eigenvalues and corresponding wave functions are calculated for the if2, HCl, CO and LiH 
diatomic molecules and AIM results are compared with the findings of the SUSY j^, the 



hypervirial perturbation method (HV) 14] , the Nikiforov-Uvarov method (NU) jll| and the 
shifted and modified shifted 1/N expansion methods [8|, |l2] as well as with the variational 
method |13]. Finally, section |TV| is devoted to the summary and conclusion. 



II. BASIC EQUATIONS OF THE ASYMPTOTIC ITERATION METHOD (AIM) 



We briefly outline the asymptotic iteration method here and the details can be found in 



references 



Q 



16l |. The asymptotic iteration method is proposed to solve the second-order 



differential equations of the form 



y" = Xoix)y' + so{x)y 



(1) 



where Ao(a;) 7^ and So(x), Ao(x) are in Coo(a,b). The variables, So(x) and Ao(x), are 
sufficiently differentiable. The differential equation ^ has a general solution [15| 



y{x) = exp 



a{x )dx 



>i -. ^ , // 



C2 + Ci / exp I {Xq{x ) + 2a(x ))dx dx 



(2) 



if /c > 0, for sufficiently large fc, we obtain the a{x) values from 

Sk{x) Sk-l{x) 



Xk{x) Xk-i{x) 



a{x), k = 1,2, 3, 



(3) 



where 



Afc(x) = X'k_i{x) + Sk-i{x) + Xo{x)Xk-i{x) 

Sk{x) = s'^_^^{x) + so{x)Xk-iix), k = 1,2,3,. 



(4) 



The energy eigenvalues are obtained from the quantization condition. The quantization 
condition of the method together with equation Q can also be written as follows 



hix) = Xkix)sk-i{x) - Xk-iix)skix) = 0, 



K — L ^ Zi^ O J . . . 



(5) 



For a given potential such as the rotating Morse one, the radial Schrodinger equation is 
converted to the form of equation (^. Then, so(a;) and Ao(x) are determined and Sfc(x) and 
Afc(x) parameters are calculated. The energy eigenvalues are determined by the quantization 
condition given by equation (0). However, the wave functions are determined by using the 
following wave function generator 



yn{x) = C2exp{- I ^ dx) 

Xk[x') 



(6) 



III. CALCULATION OF THE ENERGY EIGENVALUES AND EIGENFUNC- 
TIONS 

The motion of a particle with the reduced mass /i is described by the following Schrodinger 
equation: 



-h^ f d^ 2d 1 



2/i \dr'^ r dr r^ 



^^mV m) ^ sin" 6 86^ 



V^(r)Wn£m(r,^, 



E^r,im{r, 9, 
(7) 



Ur,Ar)Yf 



niV' Iim\ 



we obtain the 



The terms in the square brackets with the overall minus sign are the dimensionless angular 
momentum squared operator, L^. Defining \E'„^m(^, 6*, 
radial part of the Schrodinger equation: 



flj.2 r dr J "' h? 



Vir) + 



2/ir2 



, , 2uE 
Uni{r) + -^Uni[r) = 



It is sometimes convenient to define Uniir) and the effective potential as follows: 



Uniir 



'^'W, K« = nr) + '<'^'"' 



r 



2fir^ 



Since 



V dr"^ r dr J r r dr"^ 



The radial Schrodinger equation given by equation (jH)) follows that 

d'^Rniir) , 2/i 



dr 



, +^[E- Veff] Rniir) = 



(9) 



(10) 



(11) 



Instead of solving the partial differential equation ((Tj) in three variables r, 6 and 0, we now 
solve a differential equation involving only the variable r, but dependent on the angular 
momentum parameter i, which makes the solution of this equation difficult for £ 7^ or 
sometimes impossible within a given potential. 

The Morse potential we examine in this paper is defined as 



VMorseir) = D (e'^"- - 2e-"^) 



(12) 



with X = {r — re) /re and a = are- Here, D and a denote the dissociation energy and Morse 
parameter, respectively, re is the equilibrium distance (bound length) between nuclei and 
a is a parameter to control the width of the potential well. For the H2 diatomic molecule. 



the effective potential, which is the sum of the centrifugal and Morse potentials, is shown 
in figure ^ for various values of the orbital angular momentum. The superposition of the 
attractive and repulsive potentials results in the formation of a potential pocket, whose 
the width and depth depend on the orbital angular momentum quantum number for a 
given molecular potential. The potential pocket becomes shallower as the orbital angular 
momentum quantum number £ increases, which also indicates that the number of states 
supported by the potential decreases. This pocket is also very important for the scattering 
case due to the interference of the barrier and internal waves, which creates the oscillatory 
structure in the cross-section. The effect of this pocket can be understood in terms of the 
interference between the internal and barrier waves that corresponds to a decomposition of 
the scatteHng amplitude iMo two co,„po.e,.s. the ...e. a.d extent, wave. HSQ. 
The effective potential together with the Morse potential for £ 7^ can be written as. 



f(f 4. ])h'^ 
Kffir) = Veir) + VMorseir) = ^ ^^J + D [f^^^ - 26"°^) (13) 

It is known that the Schrodinger equation cannot be solved exactly for this potential for £ 7^ 
by using the standard methods such as SUSY and NU. As it is seen from equation ()13|) . the 
effective potential is a combination of the exponential and inverse square potentials, which 
cannot be solved analytically. Therefore, an approximation has to be made: The most widely 
used and convenient one is the Pekeris approximation. This approximation is based on the 
expansion of the centrifugal barrier in a series of exponentials depending on the internuclear 
distance, keeping terms up to second order, so that the effective ^-dependent potential keeps 
the same form as the potential with ^=0 [2]. It should be pointed out, however, that this 
approximation is valid only for low vibrational energy states. In the Pekeris approximation, 
by change of the coordinates x = {r — rf.)/rf., the centrifugal potential is expanded in a series 
around x = 

''''^) - '-^jrh^ - ^'(^ - ^^ + '^' ' '^' + •■■) (") 

where 7 = \^ J — . Taking up to the second order degrees in this series and writing them in 
terms of exponentials, we get 

Veix) = 7 (co + cie-°^ + cae-^"-) (15) 



In order to determine the constants Cq, Ci and C2, we also expand this potential in a series 
of X 

Vi{x) = 7 (^co + ci + C2 - (ci + 2c2)ax + (y + 2c2)«V) . . .) (16) 

Comparing equal powers of equations (fT^ and (fT^ . we obtain the constants cq, ci and C2 

as, 

. 3 3 4 6 13 ,^^, 

Co = 1 h ^, Ci = -, C2 = h ^ (17) 

Now, the effective potential with Pekeris approximation becomes, 

V,ff{x) = 7(co + cie"'^^ + C2e-2°^) + D (e^^"^ - 2e-"^) (18) 

Instead of solving the radial Schrodinger equation for the effective potential given by equa- 
tion (jl3|) . we solve the radial Schrodinger equation for the new effective potential given by 
equation ()18p obtained by using the Pekeris approximation. Inserting this effective potential 
equation (fTHj) into equation (fTT|) and using the following ansatzs 

-e'=^-^{E-jCo), /?,2 = ^(2D-7Ci), P', = ^-^{jc, + D) (19) 

The radial Schrodinger equation takes the following form: 

^^^^ + i-e' + Ple-^ - (3ie-'-^) R^,{x) = (20) 

If we rewrite equation (j^n|) by using a new variable of the form y = e~"^, we obtain 



(fRntiy) 1 dRni{y) 



a^ y2 Q,2 y Q,2 



Rni{y) = (21) 



dy"^ y dy 

In order to solve this equation with AIM for £ 7^ 0, we should transform this equation to 
the form of equation (^. Therefore, the reasonable physical wave function we propose is as 
follows 

Rni{y) = y^-e-^yfni{y) (22) 

If we insert this wave function into the equation (J2H) . we have the second-order homogeneous 
linear differential equations in the following form 

d'^fneiy) [ 2(32ay - 2ea - a^\ dfni{y) , f2e(32 + ap2-pr 



,2V 2 I ^ • \ 2 / ^«^(^) (23) 

dy \ yoi J dy \ ya'^ J 

which is now amenable to an AIM solution. By comparing this equation with equation (P), 
we can write the Ao(i/) and So(i/) values and by means of equation (jl}, we may calculate 



\k{y) and Sk{y)- This gives (the subscripts are omitted): 
^ ,2P2y-2e-a 

So 



ay 
2e(52 + al32 - Pi 



a^y 



Ai = - 

si = 2 
. . . etc 



-3P2ay + 6ae + 2a^ -6yeP2- ypi + 4 p2 y + 4 e 

a'^y'^ 
{2e P2 + a P2 - Pi^) {-a + P2y - e) 
a^y"^ 



(24) 



Combining these results with the quantization condition given by equation ^ yields 

1 ap2 - Pi 



Sq _ Si 

Aq Ai 

Ai A2 

S2 _ S3 

A2 A3 

. . . etc 



£0 



£1 



£2 



2 P2 

1 3a/?2 - Pi 

2 P2 

1 5a/?2 - Pi 

2 ^ 



(25) 



When the above expressions are generalized, the eigenvalues turn out as 

Pi - {2n + l)a/?2 



^n£ 



2P2 



n 



0,1,2,3, 



Using equation (fT9|l . we obtain the energy eigenvalues E^ 



En 



h^ 



2ftrl 



PI 
2P2 



[n + -)a 



+ 7C0 



(26) 



(27) 



As it is seen that the energy eigenvalue equation is easily obtained by using AIM. This is 
the advantage of the AIM that it gives the eigenvalues directly by transforming the radial 
Schrodinger equation into a form of y" =Xo{r)y' + SQ{r)y. In order to test the accuracy of 
equation ^I7\i. we calculate the energy eigenvalues of the H2, HCl, CO and LiH diatomic 
molecules. The AIM results are compared with those obtained by SUSY method [^ using 
original Pekeris approximation, the hypervirial perturbation method (HV) [14 1, the shifted 
1/N and modified shifted 1/N expansion methods [8y for the H2 diatomic molecule in Table 
HI In Table m we show the same comparison for the HCl diatomic molecule. Furthermore, 
the AIM results are compared with those obtained by NU method ll|, shifted 1/N and 
modified shifted 1/N expansion methods J8| for the CO and LiH diatomic molecules in 



Tables IIIII and IIVI respectively. As it can be seen from the results presented in these tables 
that the AIM results are in good agreement with the findings of the other methods. 

After we find the energy eigenvalues, the following wave function generator can be used 

to find fn{y) functions by using AIM 

/v 
'fH) (28) 

where n represents the radial quantum number and k shows the iteration number. Below, 
the first few f{y) functions can be seen 

/o(y) = 1 (29) 

/^(^)-(^«/^-^?)(^-^;^^T§t7Ty) (30) 

Hy) = i0t - ^»Pm - 3a/3.) I 1 - -^^^t-7, + ....^-5... 'l'l^-5.,. ^ I (31) 



h{y) = {-Aap2 + (3i){0t " 5a/32)(/9i - 6a/?2) 1 



o2^/p2 c„ o \/o2 n^o ^ I -, ^^^V 



/320 



"(^^^^^ + 1) 



+ „2(^?zI^ + i)(^?zI^ + 2)"a3(^?zI^ + l)(/!^^ + 2)(^^^ ^ ^ 

^ a/32 ' ^ 0-P2 ' ^ ap2 ' ^ o-Pi ' ^ cip2 ' , 

. . . etc 
It can be understood from the results given above that we can write the general formula 
for fn{y) as follows. 



fn{y) = i-ir i IT (/5i -ik+ l)«/52) ] iF,{-n, ^ + 1; ^) (33) 

\k=n / 

Thus, we can write the total radial wave function as below, 

Rni = (-1)" iUiP! -{k + l)aP,)] y'^e-^y.F.i-n, ^ + i- ^) (34) 

\k=n / 

When the hypergeometric function is written in terms of the Laguerre polynomials, we get 

Rn. = Ny'-^e-'-^yi^ (^,) (35) 

Where A^ is the normalization constant and can be obtained from 

dy = 1 as below. 



00 ^ „o r 2£„ _ n2 



N^jy'^e-'^y LT{^y) 







C+i 



-M^rv^ 



where ^ = -^ 



IV. CONCLUSION 

We have shown an alternative method to obtain the energy eigenvalues and correspond- 
ing eigenfunctions of the rotating Morse potential using Pekeris approximation within the 
framework of the asymptotic iteration method. The main results of this paper are the energy 
eigenvalues and eigenfunctions, which are given by equations i\27\\ and (jH^j) respectively. The 
energy eigenvalues are obtained for the H2, HCl, CO and LiH diatomic molecules. Our 
AIM results are compared with the findings of the other methods such as the SUSY J2[, 
the hypervirial perturbation method (HV) 'l4] , the Nikiforov-Uvarov method (NU) 11[ and 



12| as well as the variational 



the shifted and modified shifted 1/N expansion methods 
method 13| in Tables HI HIl Hill and IIVI The advantage of the asymptotic iteration method 
is that it gives the eigenvalues directly by transforming the radial Schrodinger equation into 
a form of y" =\o{r)y' + So(r)y. The wave functions are easily constructed by iterating the 
values of So{r) and Xo{r). The method presented in this study is a systematic one and it 
is very efficient and practical. It is worth extending this method to the solution of other 
interaction problems. 
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AIM SUSY HV Variational Modified Shifted 1/N Shifted 1/N 
n i results results results results expansion results expansion results 

-4.47601 -4.47601 -4.47601 

5 -4.25880 -4.25880 -4.25901 

10 -3.72193 -3.72193 -3.72473 

5 -2.22052 -2.22051 -2.22051 
5 -2.04355 -2.04353 -2.05285 
10 -1.60391 -1.60389 -1.65265 

7 -1.53744 -1.53743 -1.53743 
5 -1.37656 -1.37654 -1.39263 
10 -0.97581 -0.97578 -1.05265 

TABLE I: For the H2 diatomic molecule, the comparison of the energy eigenvalues (in eV) ob- 
tained by using AIM with other methods for different values of n and i. Potential parameters are 
D = 4.7446eF, a = IM25{A^)~^, r^ = 0.7416^°, he = 1973.29ey^° and n = 0.50391amu. 



■^.^t DO 


-^.^lUV 


-^.^l^i) 


4.2563 


-4.2590 


-4.2590 


3.7187 


-3.7247 


-3.7247 


- 


-2.2205 


-2.2038 


- 


-2.0530 


-2.0525 


- 


-1.6535 


-1.6526 


- 


-1.5374 


-1.5168 


- 


-1.3932 


-1.3887 


. 


-1.0552 


-1.0499 
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AIM Variational Modified Shifted 1/N Sliifted 1/N 
n I results results expansion results expansion results 



-4.4356 


-4.4360 


-4.4355 


-4.4352 


5 -4.3968 


-4.3971 


-4.3968 


-4.3967 


10 -4.2941 


-4.2940 


-4.2940 


-4.2939 


5 -2.8051 


- 


-2.8046 


-2.7727 


5 -2.7721 


- 


-2.7718 


-2.7508 


10 -2.6847 


- 


-2.6850 


-2.6712 


7 -2.2570 


- 


-2.2565 


-2.2002 


5 -2.2263 


- 


-2.2262 


-2.1874 


10 -2.1451 


- 


-2.1461 


-2.1194 



TABLE II: For the HCl diatomic molecule, the comparison of the energy eigenvalues (in eV) 
obtained by using AIM with other methods for different values of n and L Potential parameters 
are D = 37255cm-i, a = 1.8677(^°)-\ r^ = 1.2746^°, he = 1973.29eFylO and /i = 0.9801045amu. 
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AIM NU Variational Modified Shifted 1/N Sliifted 1/N 

n i results results results expansion results expansion results 









-11.0915 


-11.091 


-11.093 


-11.092 


-11.091 




5 


-11.0844 


-11.084 


-11.085 


-11.084 


-11.084 




10 


-11.0653 


-11.065 


-11.066 


-11.065 


-11.065 


5 





-9.7952 


-9.795 


- 


-9.795 


-9.788 




5 


-9.7883 


-9.788 


- 


-9.788 


-9.782 




10 


-9.7701 


-9.769 


- 


-9.770 


-9.765 


7 





-9.2992 


-9.299 


- 


-9.299 


-9.286 




5 


-9.2925 


-9.292 


- 


-9.292 


-9.281 




10 


-9.2745 


-9.274 


- 


-9.274 


-9.265 



TABLE III: For the CO diatomic molecule, the comparison of the energy eigenvalues (in eV) 
obtained by using AIM with other methods for different values of n and i. Potential parameters 
are D = 90540cm-i, a = 2.2994(^0)"!, re = 1.1283^°, he = 1973.29eFylO and ^ = 6.8606719amu. 
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AIM NU Variational Modified Shifted 1/N Shifted 1/N 
n i. results results results expansion results expansion results 









-2.4289 


-2.4287 


-2.4291 


-2.4280 


-2.4278 




5 


-2.4013 


-2.4012 


-2.4014 


-2.4000 


-2.3999 




10 


-2.3288 


-2.3287 


-2.3287 


-2.3261 


-2.3261 


5 





-1.6477 


-1.6476 


- 


-1.6402 


-1.6242 




5 


-1.6238 


-1.6236 


- 


-1.6160 


-1.6074 




10 


-1.5607 


-1.5606 


- 


-1.5525 


-1.5479 


7 





-1.3776 


-1.3774 


- 


-1.3682 


-1.3424 




5 


-1.3550 


-1.3549 


- 


-1.3456 


-1.3309 




10 


-1.2958 


-1.2957 


- 


-1.2865 


-1.2781 



TABLE IV: For the LiH diatomic molecule, the comparison of the energy eigenvalues (in eV) 
obtained by using AIM with other methods for different values of n and L Potential parameters 
are D = 20287cm-i, a = 1.1280(^0)-!, r^ = 1.5956A°, he = 1973.29eFA0 and ^ = 0.8801221amu. 
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FIG. 1: The shape of the rotating Morse potential for H2 diatomic molecule is plotted against the 
separation r for different orbital angular momentum quantum numbers. 



15 



